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ERRATA 

NRL  Memorandum  Report  3760 

p.  1,  bottom  line,  inside  the  parentheses,  should  read 
(ls2s  - lsjp,  ls2p  - ls^s,  ls2p  - lsjd) 

p.  2,  line  4-,  should  read  2p-3d  not  2p-3p 

p.  2,  second  paragraph,  beginning  of  line  7>  should  read, 

f N _ f N+l 
s s 

p.  14- , last  sentence  in  figure  caption  (fig.  3)  should  read: 

"Results  for  a strict  blackbody  pumping  source  give 
an  order  of  magnitude  higher  gain  at  these  densities" 
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We  show  by  means  of  a new  analytic  formalism,  which  is  developed  for  describing  the  average 
effect  of  photoexcitation  processes  in  optically  thick  media,  that  gain  can  be  achieved  in  the  AS 

XII  3-2  transitions  at  44A  at  densities  well  below,  and  over  spatial  extents  will  beyond,  those  here- 
tofore proposed.  Photoexcitation  of  the  n * 3 level  in  AS  XII  by  the  ls^  1S  — ls2p  *P  transition 
in  Si  XIII  is  used  to  generate  the  population  inversion  between  the  n - 2 and  n - 3 states  in  AS 

XII.  The  gain  is  enhanced  in  magnitude  and  spatial  extent  by  photoexcitation  of  the  ls2p  *P 
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20.  Abstract  (Continued) 

and  ls3p  1P  states  of  the  Si  and  A2  plasmas  respectively.  An  estimate  of  the  inherent  experimental 
limitations  of  this  inversion  technique  is  also  made  by  comparing  the  calculated  gains  with  upper 
limits  that  are  obtained  by  pumping  with  a strict  blackbody  photon  flux  in  the  pumping  line. 
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I . ^Introduction 

"^Lt  is  a well  established  phenomenon  of  radiative  transfer  in  stellar- 
like  plasmas  that  as  the  optical  depth  of  a line  increases,  the  mean 
radiative  energy  flux  in  the  line  correspondingly  increases  to  approach 
the  Planck  blackbody  value.  Recently,  the  possibility  of  achieving 
stimulated  emission  in  the  x-ray  or  soft  x-ray  region  by  utilizing  the 
interaction  of  strong  lines  from  different  elements  or  neighboring  ioni- 
zation stages  has  been  considered  in  the  context  of  resonance  line  pump- 
ing  in  high  density,  albeit  optically  thin,  plasmas.  * -We-  demonstrate 
-hero- that  proper  utilization  of  photon  trapping  effects  in  the  pumped 
and  pumping  plasmas  can  yield  comparable  gains  at  soft  x-ray  frequencies 


to  those  proposed  -is 


however,  at  plasma  densities  that  are 


significantly  below  and  over  spatial  extents  that  are  significantly 
larger  than  those  that  would  be  required  were  one  to  be  restricted  to 
optically  thin  plasmas  in  which  an  exponentially  decaying  inversion 
density  is  generated.  The  demonstration  is  made  by  treating  the  physi- 
cal situation  of  a non-collisionally  dominated  aluminum  plasma  with  the 
aid  of  a new  global  rate  equation  formalism  developed  to  include  photo- 
excitation processes,  on  the  average,  into  ionization  equilibrium  calcu- 
lations involving  optically  thick  laboratory  plasmas.  -In  particular, 
the  pumping  of  the  kl  XII  ls^  - ls3p  singlet  transition  (6^655A)  by  the 
Si  XIII  Is  - ls2p  principal  resonance  line  at  6.650!  is  considered  *or 
a range  of  plasma  densities,  sizes,  and  mass  ratios  with  a specific 
analysis  of  photon  trapping,  pumping,  and  escape  within  both  the  silicon 
and  the  aluminum  plasmas.  Of  the  three  A £ XII  singlet  transitions  in 

which  gain  may  be  anticipated  (ls2s  - ls3p  - ls2p,  - ls3s,  ls2p  - lsjd), 
Note:  Manuscript  submitted  March  23,  1978 
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the  2p-3d  is  the  strongest  transition;  hence  it  will  have  the  largest 
gain  under  the  assumption  we  make  that  the  s,  p,  and  d states  of  levels 
2 and  3 are  statistically  populated.  In  this  case,  the  gain  in  the 
2s-3p  line  will  be  about  20^  of  the  gain  in  the  2p-3p  line,  while  for 
the  2p-3s  line,  it  is  about  3$ 


II . Pumping  Model 

Consider  first  the  origin  and  the  random  walk  of  pump  photons 

through  an  optically  thick  silicon  plasma.  If  the  single-flight  photon 

escape  probability  is  P , and  the  probability  for  collisional  quenching 

per  re-absorption  is  P then  the  probability  that,  in  scattering,  the 

photon  neither  escapes  nor  is  quenched  is  fg,  given  by  (1-P^)  (1-P^). 

The  probability  P^  that  the  photon  executes  exactly  N scatterings  is 
N N-l 

fg  - f and  thus  Ng,  the  mean  number  of  scatterings  in  the  silicon 
plasma,  is  Z)  N P^  = fg/(l-fg).  Each  scattering  represents  a re-exci- 
tation  of  the  Is2  - ls2p  transition;  hence  the  collisional  excitation 


rate  of  the  ls2p  state  is  effectively  enhanced  by  the  factor  (1  + N ) 
in  accounting  for  the  net  effect  of  both  radiative  and  the  original 
collisional  excitation. 

In  general,  this  technique's  accuracy  in  the  mean  depends  only  on 

the  accuracy  with  which  P and  P can  be  calculated  for  each  optically 

e q 

thick  transition.  The  values  of  P , which  are  given  by  c 

3 A 


ij 


+ C 


where  is  the  Einstein  A for  the  transition  and  C is  the  total 
collisional  transition  rate  out  of  the  upper  level,  were  obtained  from 
distorted-wave  calculations  of  the  collision  strengths^.  Their  values 
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are  listed,  along  with  the  Si  - five  level  structure  being  studied, 

in  Figure  1.  For  collisional  ionization  out  of  the  excited  states  the 

k 

following  expression  was  employed 


N 10 

C = e 

coll.  ion. 


-5 


(kTe/X)7 
(5  + kT  /x) 

e 


exp  C - X ! kT  ) sec 
e 


-1 


(1) 


where  X and  kT^  are  in  eV,  x being  the  energetic  gap  between  the  level 
and  the  continuum.  Note  that  to  account  for  both  spontaneous  radiative 
and  collisional  transitions  out  of  a particular  level,  the  spontaneous 
rate  can  be  divided  by  (1-P^),  the  ratio  of  radiative  to  total  transi- 
tions out  of  the  level.  It  should  be  noted  that  the  formalism  described 
here  neither  requires  nor  includes  a model  of  the  complete  ionization 
structure  of  the  aluminum  or  silicon  plasmas.  However,  we  have  made 
reasonable  assumptions  about  the  densities  of  the  active  ions,  A l XII 
and  Si  XIII,  and  the  electron  densities,  and  perform  the  computation 
using  the  rates  of  processes  occuring  with  these  ions.  For  the  initial 
clarity  of  working  with  a simple  set  of  rate  equations  we  have  not 
included  recombination  processes  involving  A l XIII  and  Si  XIV.  This 
will  impact  estimates  of  the  gain  achievable  in  the  kl  XII  3-2 
transitions.  Collisional  recombination  processes  tend  to  increase  the 
populations  of  the  upper  levels  and  thus  would  directly  enhance  the  gain 
of  the  3-2  transitions.  We  calculate  that  the  dielectronic  recombination 
coefficient  into  the  n = 2 levels  of  A l XII  is  12$  greater  than  that  into 
the  n = 3 levels.  However,  as  shown  in  figure  1,  the  electron  colli- 
sional excitation  rate  into  the  n = 2 levels  is  7 times  as  great  as 
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that  into  the  n = 3 levels  from  level  1 at  the  8C0  eV  t 'mperature  con- 
sidered here.  Since  we  show  below  that  substantial  gain  is  achievable 
with  even  this  large  collisional  population  ratio  ^precisely  because 
photo-excitation  is  dominant],  the  techniques  used  here  therefore  give 
first  order  results  by  not  including  these  processes. 

Finally,  the  single-flight  escape  probabilities  were  calculated 
from  the  following  formula  due  to  Athay^, 


P 

e 


(2) 


where  a = A^AffAv^,  is  the  line  broadening  parameter  appropriate  for 
a Voigt  profile  that  is  naturally  broadened  in  the  line  wings  and 

= the  plasma  optical  depth  at  the  line  center,  along  a plasma  radius, 
corrected  for  stimulated  emission.  This  formula  has  been  found  to  agree 
within  20$,  for  values  of  a less  than  C.C5,  with  detailed  multi-frequency 
transport  calculations  . 

For  the  calculations  presented  here  we  assumed  a cylindrical  geome- 
try with  the  cylinder  Z-axis  much  longer  than  the  radius.  The  outer 
portion  of  the  cylinder  is  silicon,  with  a temperature  of  1 Kev.  The 
core  is  aluminum,  assumed  to  be  at  C.3  Kev.  These  teraperatues  are  shown 
by  coronal  models  to  be  those  at  which  the  kl  XII  and  Si  XIII  ion 
species  constitute  about  half  of  all  the  ions  in  the  plasma.1 

. Direct  collisional  processes  involving  any  mixed  Si  XIII 

and  kl  XII  at  the  boundaries,  such  as  charge  exchange  re- 
actions, are  not  considered  since  the  diffusion  distance  is 
only  2 microns  per  nanosecond  at  the  lowest  densities  to  be 
considered  (10"  ions/cm3).  Unwanted  mixing  of  the  two  sub- 
stances could  in  any  event  be  suppressed  by  operating  them 
with  a neutral  filler  such  as  polyethelene. 
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Consider  next  the  escape  of  photons  from  the  silicon  region.  Some 

photons  will  cross  the  outer  surface  of  area  Ag  and  others  will  escape 

through  the  inner  silicon  area  (A&)  into  the  aluminum.  The  fraction 

that  enter  the  aluminum  is  thus  obtained  by  the  ratio  of  available  escape 

surfaces  A /(A  + A ) in  first  approximation.  Of  those  entering  the 

as  a 

aluminum  region,  only  a fraction  0 will  be  trapped.  If  the  aluminum 

region  is  many  absorption  lengths  thick  at  line  center  nearly  all  the 

entering  photons  will  be  intercepted  in  the  aluminum  region  (i.e.  3 a 1). 

This  is  the  physical  situation  we  have  assumed  in  our  calculations. 

This  condition  will  be  more  easily  realizable  if  the  aluminum  and  silicon 

plasmas  are  counter-streaming  to  being  the  pumped  and  pumping  transition 

wave  lengths  into  coincidence  (at  a counter-streaming  velocity  of 

7 x 1C7  -^7  ).  In  general,  however,  the  total  probability  € that  a 

silicon  photon,  which  ultimately  escapes  the  silicon  region  with  proba- 
co  n p 

bility  P = P f 1 = i — —7 — — — r-  enters,  becomes  trapped 

u e N=0  3 (1"V 

in,  and  pumps  the  A 2 XII  Is2  - ls3p  transition  is  given  by  P^ 


By  taking  into  acoount  the  above  considerations,  and  by  assuming 
statistical  population  of  the  A 2 XII  2-  and  3“  sublevels,  one  can  aver- 
age over  the  geometrical  details  of  a radiation  transport  calculation 
and  write  a simple  set  of  rate  equations  for  the  A l X7I  steady-state 
level  populations  in  terms  of  quenching  probabilities  and  mean  numbers 
of  photon  scatterings.  Numbers  1-5  will  refer  to  the  A;  XII  levels  1, 
2,  and  3»  and  Si  XIII  levels  1 and  2 respectively.  Then,  the  main 
features  of  the  level  2 and  level  3 photo-excitation  dynamics,  based  on 
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the  level  system  shown  in  Figure  1 (AjJ  XIII  and  Si  XIV  ground  states  also 
included),  can  be  described  by  the  equilibrium  equations: 


N1  °12  (1  + = 


sl2 


1 - P, 


(3) 


Q2 


»1  °13  (1  + + \ ^5  NS1J  - ■ 


(4) 


03 


where  the  IT's  refer  to  total  numbers  in  the  plasma  rather  than  densities. 
By  solving  Eq.  (3)  and  (4)  for  the  population  ratio  N^'N^,  one  obtains 


i _ A21(1-V  [N1  C13  (1  + Ns13}  + \ C45  gNsl3 
N2  A,,  N,  C10  (1+N10)  (1  - Pno) 


(5) 


31 


12 


sl2 


Q2 


N e 

An  inversion  exists  when  3/N2  > 63/g2  = In  these  equations,  the 

C's  are  collisional  excitation  rates  for  the  indicated  transitions,  and 
the  A's  are  spontaneous  radiative  rates.  Ng  refers  to  the  calculated 
mean  number  of  scatterings  in  the  Xl  XII  region  for  the  indicated  tran- 
sition. Equation  (5)  clearly  exhibits  the  influence  of  the  trapped  line 
radiation  on  the  pumping  process;  population  inversion  between  levels  2 
and  3 is  enhanced  by  scatterings  to  level  3 and,  naturally,  reduced  by 
scatterings  to  level  2 which  sustain  population  in  the  lower  level  of  the 
lasing  transition.  As  expected,  the  more  Si  XIII  ground  states  that 
can  be  excited,  the  larger  is  the  attainable  inversion.  Since  we  do 
not  attempt  in  this  treatment  to  account  for  radiative  depletion  of 
level  3 once  lasing  begins,  Eq.  (5)  will  describe  only  the  linear 
amplifier  behavior  of  the  aluminum  plasma. 
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III.  Numerical  Results  for  Gain 


From  Eq.  (5),  one  can  calculate  the  volume  averaged  steady-state 
gain  coefficient  which  will  be  established  following  diffusion  of  the 
silicon  radiation  into  the  aluminum  region.  In  Figs.  2 and  3 the  gain 
coefficient  is  plotted  for  the  strongest  of  the  AH  XII  2-3  singlet 
transitions  (the  2p-3d) . Since  roughly  half  the  ions  are  in  the  A£  XII 
or  Si  XIII  states  the  electron  density  was  taken  equal  to  roughly  20 
times  the  ground  state  densities  of  these  species. 

A meaningful  upper  limit  to  the  magnitude  of  the  inversion  which 
may  be  produced  by  this  mechanism  may  be  obtained  by  assuming  that  at 
most,  the  AH  XII  region  is  optically  thick  to  a broad  band  flux  of 
radiation  that  is  10  Doppler  widths  on  either  side  of  line  center.  In 
each  of  the  Figures  2 and  3 we  have  indicated  the  gain  which  would  be 
produced  by  a strict  blackbody  flux  of  pumping  photons  at  1 keV,  20 
Doppler  widths  wide,  impinging  on  the  A H XII  region.  For  the  laboratory 
plasma  conditions  considered  in  this  paper,  it  is  always  much  greater 
than  the  gain  calculated  from  Eq.  (5).  Of  course  if  the  A;,  XII  region 
is  not  thick  for  quite  20  Doppler  widths,  this  "upper  limit"  will  be  in 
fact  smaller  since  many  of  the  blackbody  photons  would  simply  pass 
through  the  AH  XII  region. 

Fig.  2 illustrates  that  the  gain  in  AH  XII  is  dramatically  increased 
by  increasing  the  Si  XIII  density;  this  ultimately  provides  much  more 
pumping  radiation  by  increasing  the  total  number  of  collisionally  created 
Si  XIII  photons.  The  ratio  of  total  Si  XIII  to  AH  XII  ground-state  ions 
was  kept  at  5^  in  this  particular  calcu  ition.  Figure  3 illustrates  a 
gain  saturation  effect.  For  a Si  XIII  to  A l XII  ground  state  ion  ratio 
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of  54,  the  gain  is  maximized  at  an  A l XII  number  density  of  about 
19  -5 

1 x 10  cm  . Above  this  density,  collisional  quenching  of  the 

Af  XII  1-3  singlet  transitions  reduces  the  gain  considerably;  we  also 

expect  significant  gain  gradients  to  develop  at  higher  densities  as 

pumping  photons  from  the  silicon  region  find  it  increasingly  difficult 

to  penetrate  to  and  pump  the  core  of  the  A2  XII  region.  Below  this 

density  the  reduced  density  of  lasing  ions  is  responsible  for  the 

reduced  gain.  We  also  note  that  significant  gain  is  achieved  at  alumi- 

13  “3 

num  ion  densities  as  low  as  10  cm  , which  is  more  than  four  orders 
of  magnitude  lower  than  the  densities  suggested  in  reference  1. 

Finally,  the  quenching  effect  of  radiation  trapping  in  the  A i XII 
1-2  transition  can  also  be  seen  in  Fig.  2.  The  steady  state  gain  co- 
efficient for  a fixed  density  and  ion  number  ratio  is  plotted  against 
the  radius  of  the  aluminum  region.  The  gain  coefficient  is  monotoni- 
cally  increased  by  increasing  the  aluminum  region's  size;  since  S has 

been  taken  equal  to  one,  this  is  purely  an  effect  of  the  increased 

2 

optical  depth  of  the  A£  XII  region  in  the  A£  XII  Is  - ls2p  principal 
resonance  line. 


IV.  Summary 

In  plasmas  whose  excitation  dynamics  are  dominated  by  radiation 
transport  the  degree  of  photo  excitation  attainable  can  be  calculated 
in  the  mean  from  locally  definable  collisional  quenching  and  nonlocal ly 
definable  photon  escape  probabilities.  The  escape  probabilities  are 
functions  of  both  the  plasma  geometry  and  the  absorption  and  emission 
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profiles.  However,  in  simple  geometries  and  for  Doppler  or  Voigt  pro- 
files, simple  expressions  for  the  escape  probabilities  as  a function  of 
the  mean  optical  depth  at  line  center  can  be  derived  from  detailed  two- 
level  radiation  transport  calculations.  In  this  paper,  we  showed  how 
these  probabilities  could  then  be  utilized  in  a global  rate  equation  to 
determine,  in  steady  state,  the  mean  population  inversion  and  gain  co- 
efficient achievable  in  a two  component  Si  - A L plasma  for  different 
sizes,  densities  and  mass  ratios  in  a cylindrical  configuration  in  which 
the  aluminum  was  surrounded  by  the  silicon  plasma. 

A technique  for  producing  at  least  some  of  these  configurations  is 
that  of  laser  heating.  If  the  outer  silicon  region  is  to  be  heated  to 
1 Kev,  and  the  inner  aluminum  region  to  0.8  Kev,  an  approximate  expres- 
sion for  the  heating  time  of  the  aluminum  core  by  electron  thermal  con- 
duction alone  may  be  obtained  by  estimating  the  temperature  gradient 
during  heating  to  be  1 Kev/R^,^  where  Is  the  outer  radius  of  the 

cylinder.  (We  assume  the  aluminum  to  be  preionized  by  an  electron  beam, 

for  example.)  Using  the  expressions  for  electron  conduction  given  by 

7 

Braginskii  one  finds  that 

'heat  “ * 10'27  Rtot  “ <6) 

where  t,  _ is  the  time  in  seconds  needed  to  heat  the  center  to  0.8  Kev, 
R(A2)  is  the  radius  of  the  aluminum  region  in  centimeters,  N is  the 
average  ion  density  (taken  to  be  0.05  times  the  electron  density)  for 
ki  XII  and  Si  XIII.  Eq.  (6)  indicates  that  it  would  be  possible  to  heat 
the  aluminum  interior  by  thermal  conduction  before  the  silicon  plasma 


disassembles  (at  t 1 nsec)  for  characteristic  radii  of  C.l  cm  and  ion 

ig 

densities  a few  times  1C  or  less.  For  larger  radii  and  higher  den- 
sities, the  core  of  aluminum  will  have  to  be  heated  directly  by  the 
laser  or  by  another  technique  such  as  in  a pulsed-power  exploding  wire 
experiment . 
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Si  XIII  Al  XII 

SINGLETS 


TRANSITION 

A 

C (per  ion  per 
electron  per  sec) 

1 

3.9  X1013  sec1 

5.5  X10'12 

2 

2.8  X1013 

6.0  X10'12 

3 

1.4  X1012 

3.9X10  10 

4 

4.6X1011 

1.0X10'10 

5 

2.2X1011 

6 

7.7X1012 

8.5X10  13 

Fig.  1 - The  energy  levels  of  interest  in  the  AC  Xll-Si  XIII 
lasing  system  are  indicated  in  this  diagram  along  with  the 
radiative  and  collisional  transition  probabilities  used  in 
the  calculations. 
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